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ON THE GENERALIZED BUCKLEY-LEVERETT 

EQUATION 

JAN BURCZAK, RAFAEL GRANERO-BELINCHON, AND CARVING K. LULI 


Abstract. In this paper we study the generalized Buckley-Leverett 
equation with nonlocal regularizing terms. One of these regularizing 
terms is diffusive, while the other one is conservative. We prove that 
if the regularizing terms have order higher than one (combined), there 
exists a global strong solution for arbitrarily large initial data. In the 
case where the regularizing terms have combined order one, we prove 
the global existence of solution under some size restriction for the initial 
data. Moreover, in the case where the conservative regularizing term 
vanishes, regardless of the order of the diffusion and under certain hy¬ 
pothesis on the initial data, we also prove the global existence of strong 
solution and we obtain some new entropy balances. Finally, we provide 
numerics suggesting that, if the order of the diffusion is 0 < a < 1, a 
finite time blow up of the solution is possible. 


Contents 


1. Introduction 

1.1. Aim and outline 

1.2. Results 

1.3. Notation and preliminaries 

2. The entropy inequalities 

3. Proof of Proposition 1; Weak solutions 

4. Proof of Proposition 2: Decay estimates 

5. Proof of Theorem 1; Global solutions for max{a,/3} > 1 

6. Proof of Theorem 2: Global solution for max{a, /3} = 1 

7. Proof of Theorem 3: Global solution if 0 < a < 1 and fi = 0 

8. Proof of Proposition 3: Finite time singularities 

9. Numerical simulations 
References 


1 

2 

3 

7 

9 

11 

13 

14 

17 

18 
20 
21 
23 


1. Introduction 


In this paper we study the case of the Buckley-Leverett equation with 
generalized regularizing terms provided by fractional powers of the laplacian 


(-A)“/2 = A“ 

( 1 ) dtu + dx 


V? -h M(1 - tt)2 



0 < tt(x,0) = uq{x) < 1 


1 


with initial data 
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and where M > 0 is a fixed constant. Here H is either H = M or H = T. 

Let us immediately emphasize that uq{x) < 1 is not a smallness condition, 
since, in applications, u denotes a certain proportion (compare the following 
literature outline). Equation (1) is a nonlocal regularization of the classical 
Buckley-Leverett equation 


( 2 ) 


dtu + dx 


V? + M(1 



0, X G t > 0, 


The nonlinearity in equation (1) is regularized in two different ways: firstly, 
due to the diffusive term 




and secondly due to the conservative term 


Equation (2) was derived by Buckley &: Leverett in[4] and it has been well 
studied since then (see LeVeque [27] and Mikelic Sz Paoli [30]). This equation 
is used to describe a two-phase flow in a porous medium. For example, oil 
and water flow in soil or rock. In this situation u represents the saturation 
of water and M > 0 is the water-over-oil viscosity ratio. Equation (2) is a 
prototype of conservation laws with convex-concave flux functions (see for 
instance Lax [26], Glimm [15], Hong [22], Hong & Temple [23], Bayada, 
Martin Sz Vazquez [3]). Under the effect of dynamic capillarity, (2) needs to 
be modified with two regularizing terms (see Hassanizadeh &: Gray [20, 21]): 


( 3 ) 


dtu + dx 


V? + M(1 



izdxU + v'^rd'^dtu, x G H, t > 0. 


Equation (3) has been studied by many authors. For instance, Van Duijn, 
Peletier & Pop [31] derived existence conditions for solution of travelling 
wave form. Moreover, this leads to admissible shocks for (2), which violate 
the Oleinik entropy condition. In [24], Hong, Wu &: Yuan proved the global 
existence of classical solution to (3) with uq G iL^(M). Furthermore, they 
proved that the solution becomes C'°°((0,oo) x M). For equation (2), Hong, 
Wu Sz Yuan proved that if the total variation T.V.{uq) is sufficiently small, 
then there exists a solution to (2). Wang & Kao [32] studied (3) on a 
finite interval H = (0, L) and showed that the solution, ul, converges to the 
solution Uoo on the half line as L —)• oo. 


1.1. Aim and outline. The purpose of this paper is to study (1). We are 
mainly interested in global existence of solution together with their qualita¬ 
tive behaviour as well as in the finite time singularities. 

We provide details of our results in the following subsection. Subsection 
1.3 contains notation, including the definition of a weak solution and certain 
preliminaries. Section 2 provides new entropy inequalities for the fractional 
laplacian that are interesting by themselves, therefore these inequalities are 
stated for an arbitrary dimension d. Sections 2-8 contain proofs of our re¬ 
sults. Finally, in Section 9, we provide some numerical results suggesting 
the existence of finite time singularities for the cases 0 < a < 1 and /r = 0. 
These numerics suggest also that in the critical case a = 1 the solution 
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exists globally. This is in agreement with the results for the Burgers equa¬ 
tion with fractional dissipation by Kiselev, Nazarov & Shterenberg, [25] and 
Dong, Du & Li [14]. Let us remark that, when the term fiA^dtu is added 
to the equation, even for a = (5 = 0.25 there is no evidence of blow-up. 
Consequently, our numerics appear to discard a finite time blow-up scenario 
when /r > 0. 

To the best of our knowledge, all our results are new. 


1.2. Results. First, let us provide a result concerning the global existence 
of weak solutions for (1), corresponding to rough initial data, i.e., merely 
0 < uq < 1 a.e., as well as concerning new entropy balances (that are needed 
in the existence part of the result, but are interesting by themselves). 

Proposition 1. Let 0 < uq < 1, uq G L^(D) nL°°(D) be the initial data for 
(1) with u>0, 0<a<2, /i = 0 and M > 0. Then there exists a global 
weak solution such that 


u G L°°(0, oo; L^(D) n L“(L!)) n oo; 

Furthermore, if u is an L^(0, T; R^(D)) solution to (1), then the following 
entropy inequalities hold 

(4) / u{t)\og{u{t))dx + i' / / A°‘u{s)log{u{s))dxds < / mq log(uo)c? 2 :, 

Jn Jo Jn Jn 


and 

(5) 


u{t) + 


M 


1 + M 
2M^-^ 


l + M 

rt 

+ V 


log (n(t)^ -b M(1 — n(t))^) dx 

(VmI ^ 


I 

Jn 


arctan 


— 1 ] dx 


-b 


\u{t) 

f f A“u(s) log (m(s)^- b M(1 — n(s))^) 
Jo Jn i 

2{M + l)u{sf 


< 


u{sy -b M(1 — tt(s))2 
M 


dxds 


uo + 


- 2 - 


m 

M 


l + M 


log {ul -b M(1 - uo)^) dx 


1.5 


1 -b Af 


arctan ( y/M f -1 ) ) dx. 

: V \uo JJ 


Let us remark that the terms 

A"rt(s) log(u(s))dxds, 



0 Jn 


and 


A°u(s) 


0 Jn 


log (u(sf + M(1 - u(s)f) + 


2(M + l)u(sy 
u(s)^ -b AL(1 — u(s))^ 


dxds 


will provide a bound on a fractional derivative of the solutions. The 
proof of Proposition 1 will be established in Section 3. 
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Our next results concern the qualitative behaviour of smooth solutions. 
In case Q = T, we denote the average of uq by 


(«o) 



Proposition 2. Let u be the classical solution to (1) with initial data 0 < 
Uq < 1; where u > t), < a <2 and M > 0. Then, 

(1) for n = T.- 

2r(l+Q:) cos((l —Q:)7r/2) , 

II^^WIIl°°(t) < {uo) + (llnollioo - (uo))e _ 


(2) for n = M.- 


\u{t)\\L^ < 


Ikolli 


1 -j- o 


r(l+a) cos((l—a)7r/2) 


2n 


Our main results address the problem of global existence of smooth so¬ 
lutions. More precisely, we have results for three cases, depending on the 
values of the parameters a and /3. 

(1) the subcritical case: the higher space derivative is in the dissipative 
term, i.e. 1 < max{a,/I} < 2. Here we show global existence of 
smooth solutions with no restrictions on the initial data. Compare 
Theorem 1. 

(2) the critical case: the transport term exactly balances the regularizing 
terms, i.e. 1 = max{a,/3}. Here, for > 0 and /3 = 1, we prove 
global existence of smooth solutions without any size restriction on 
the initial data. In the other cases we need certain smallness con¬ 
ditions. Namely, for p. = 0, u > 0 and a = 1, we obtain global 
existence of smooth solutions for initial data satisfying a smallness 
restriction on the lower order norm this smallness restriction is 
explicit in terms of iz and M. Finally, in the case /r > 0, a = 1 and 
0 < /3 < 1, we obtain the global existence for initial data satisfy- 

1+/3 

ing a smallness condition in FI 2 . The smallness restriction is here 
slightly less explicit, but easily computable. See Theorem 2. 

(3) the supercritical case: the higher space derivative is in the transport 
term, i.e. 0 < a < 1 and /r = 0. Even here, for = T, we 
are able to prove global existence of smooth solutions for smooth, 
periodic initial data satisfying an explicit smallness restriction on 
the Lipschitz norm 

The remaining open problems are in the critical and supercritical regime. 
In particular, our results do not apply to the case where 

max{a,/?} < 1, > 0, 

and there is no large data, global results for the critical case with p = 
0,12 > 0,a = 1. In the context of the latter, let us observe that on one 
hand, there are certain new methods available for nonlinear problems with 
nonlocal critical dissipation, like the method of moduli of continuity by 
Kiselev, Nazarov &: Shterenberg [25], the hne-tuned DeGiorgi method by 
Caffarelli & Vasseur [6] or the method of the nonlinear maximum principles 
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by Constantin &: Vicol [9] (see also Constantin, Tarfulea & Vicol [10]). But 
on the other hand, our nonlinearity is more complex than the typical ones. 

Now, let us state the main theorems. First, we study the subcritical case 
max{a, /?} > 1 . 


Theorem 1. Let 0 < uq < 1, uq G be the initial data for (1) with 

M > 0. Then (1) has a global solution 

u G C([0,T], n l2(0, T; for allO<T <oo. 


Moreover, for t <T, the solution satisfies 


W{t)\\l2 + /r||w(t)||^^/2 + 2u [ ||u(s)||^„/2fis = ||ixo|li2 + fi\\uo\\%^/2, 

Jo 


provided 

• either u> 0 , l<a< 2 ,/r = 0 and s > 1 (purely dissipative regular¬ 
ization), 

• or /r > 0 , u > 0 , 1 < max{Q;,/?} < 2 and s > max{l + 

(dissipative-conservative regularization). 


For the critical case, let us define the following constants 


Definition 1. Let 7 * be a constant such that 

2j*{M + l) , 2(7*)2(7*+ M)(M + 1)2 _ 


M 


M2 


and let 7 he any fixed number such that 0 < 7 < 7 *. 

Next, let Cs he the Sobolev’s constant corresponding to the embedding 

-A L°°. 


We have 

Theorem 2. Let 0 < uq < 1, uq & M^(n), s > 1 be the initial data for (1) 
with M > 0. Then (1) has a global solution 

u{t) € C([0, T], M*(fl)) n l2(0, T; V T < 00 

that satisfies the energy balance 

Mt)\\l2 + M||u(t)||^^/2 [ h(s)||^./2fis = ||uo|li2 + /r||no||^^/2. 

Jo 

Under the following conditions: 

(i) Either: u>0, 0<a<2 and fi > 0, fi = 1 (conservative regulariza¬ 
tion with no smallness conditions on the data). 

(ii) Or: iz>0, a = l, fi = 0 and the initial data is such that 

\\uo\\l^ < 7 - 

In this case the solution satisfies the maximum principle 
||'w(i)llH'0-5 < II^^o||hO-5- 


(6) 
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(iii) Or: a = l,/r>0, 0</3<l and the initial data is such that 

(7) lko|li2 + (1 + h)\Wo\\%o.5 + /^ll'^oll^l^ < (1 + 

s 

Then, the solution satisfies the maximum principle 

ll'^(^) llj:/0.5 T /i||rt(t)|| < ll'^oll^O.S “I" A^ll'^oll ■ 1+|3 • 

H 2 U 2 


Remark 1. The fact that 


( 8 ) 


2u 

u2 + M(l-ti)2 


u^{2u — 2 M (1 — u)) 
(^2 + M(1 — ufi)^ 



< C{M), 


independently of the value of u, allows for a global result relying on a condi¬ 
tion related to M, v and //. However, we are interested in results that deal 
with every possible value of the physical parameters present in the problem. 


In our opinion, there are two reasons, at least in the case /i = 0, why 
the smallness condition ( 6 ) may be seen as a rather mild restriction. The 
first one is that the size restriction affects a lower norm, merely L°°, keeping 
the higher seminorms as large as desired. The second one is that, given M 
and v, the constant 7 * can be easily computed. For instance, if we further 
assume 7 * < 1 , the expression for 7 * is explicit: 


7 


* 


min 


— 1 + 

l + M 



In particular, 7 * = 0{n) for <C 1. 

The last case, namely where 0 < a < 1, is harder because the leading term 
in the equation is the transport term. However, under certain conditions, we 
can prove the global existence of solutions. Before we can state the relevant 
result, we need some notation. 


Definition 2. Let 7 and M be given, positive constants. Define B( 7 ) as 
follows 

, 27 (M + 1) 47 ( 7 + M)(M +1)2 472 ( 7 + M)(M +1)2 

= M +--+- Jp - 

, 273(M + 1)3 87 ^( 7 + M)2(M +1)3 

+ W + W ■ 

Next, let 7 * be a small enough constant such that 

(10) n-n < /a+“)<=”((! 

TT 

Then we have the following 


Theorem 3. Let H = T and 0 < uq < 1, uq £ with s > 2 be the 
initial data for (1), where M > 0 and > 0, 0 < a < 1, p. = 0. Fix any 
0 < 7 < 7 *, with 7 * in accordance with Definition 2. If the initial data is 
such that 

||'Mo||wl’°°(T) < 7) 

then there exists a global solution 

u G C([0,r],R*(T)) nL2(0,T;R*+"/2(Tr)) V T < 00 . 
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Furthermore, this solution satisfies the maximum principle 

||t6(t)||vKi.oo(T) < ll^^o||vKi.°°(T)) 

and the energy balance 

Mt)\\l2 + U / ||w(s)||^„/2(is = \\uo\\l2. 

Jo 

In the above theorem, we impose domain restrictions and stronger small¬ 
ness assumptions. These domain restrictions are due to the better behavior 
of the fractional laplacian in a bounded domain. The size restrictions on 
data are again on a lower order norm (Lipschitz) and with a rather explicit 
constant. 

Finally we obtain the standard finite time blow up for certain initial data 
in some Holder seminorm: 

Proposition 3. Fix a constant M > 0 and consider /r = 0, min{u, a} = 
0. Then, there exist 0 < uq < 1 £ F[‘^{Q) and T* < oo such that the 
corresponding solution, u{t), of equation (2) has a finite time singularity in 
for 0 < (5 <C 1, i.e., 

limsup ||u(t)||c'^ = oo. 
t^T* 

The proof of this result is obtained by a virial-type argument. How¬ 
ever, we remark that it can also be obtained by means of pointwise ar¬ 
guments (see Castro &: Cordoba [7] for an application of these pointwise 
arguments to prove blow up). These virial-type arguments have been used 
for several transport equations even in the case of nonlocal velocities (see 
Cordoba,Cordoba & Fontelos [12], Dong, Du &: Li [14], Li &: Rodrigo [28] 
and Li, Rodrigo & Zhang [29]). In this case, the transport term is highly 
nonlinear and this method fails in the case of viscosity u>0, 0 <q:<C1. 


1.3. Notation and preliminaries. 

1.3.1. Singular integral operators. We denote the usual Fourier transform of 
u by u. Given a function u : D ^ M, we write A“n = (—A)"/^tt for the 
fractional laplacian, i.e. 

^(0 = ier«(e)- 

This operator admits the kernel representation 


( 11 ) 


A“n(x) = ^ P.V. 

7eZ‘* 




u{x) - u{y) ^ 

X — y — 2777 ]'^+" 


if the function is periodic and 

(12) A“n(x) = C„,rfP.V. [ 


u{x) - u{y) 


\x — 7 /]'=*+“ 

if the function is flat at infinity. Notice that we have 

r(l -|- a) cos((l — 0 ) 77 / 2 ) 


dy, 


(13) 


Ca,l — 


TT 


where r(-) denotes the classical F function. 
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1.3.2. Functional spaces. We write for the usual L^-based Sobolev 

spaces with norm 

Il/lll/^ = ll/lli2+ 11/11^., II/IIh^ = IIAVIIl- 

The fractional L^-based Sobolev spaces, are 

W^’P = 1/ G LP{n^),di^^f G £ iP(^Qd X 


with norm 

where 


= wfwi. + ii/ir^..„ 


1.3.3. Entropy functionals. For a given function u > 0, we define the follow¬ 
ing entropy functionals 

(14) Fi[u]= / (u(x) log(u(x)) — u(x)-|-1) dx, 


(15) 


-F 2 N = / (1 + u(x)) log(l-|-u(x))dx. 
Jn<i 


These two entropies have an associated Fisher information: 


(16) 


M = / A“u(x) log(u(x))dx, 
Jnd 


(17) 


^2 M = / A"u(x) log(l-|-u(x))dx. 
Jnd 


The third entropy that we are using reads 

(18) F^[u]= f u{x) log [u{x)‘^ + M{1 — u{x))^) dx, 


with its Fisher information 
(19) 


h.'^uu^ 


IfH = A“«log {«" + M{1- uf) dx + 2(Af + 1 ) ^ 

1.3.4. Notation. Recall that we denote the mean of a function by 

(u) = ^ 

Let us introduce / and a as follows 

u(x, t)^ 


dx. 


( 20 ) 

and 

(21) 


f{u{x,t)) = 


u{x,ty + Af(l — u{x,t)y’ 


df 2 u{x,t)M{l - u{x,t)) 
a{u{x,t)) = 3- = 


du {u{x,ty + M{1. — u{x,t)yy 
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Finally, let us introduce the notation for the mollihers. For e > 0, we write 
jTg for the heat kernel at time t = e and define 

Je*f = Jef. 


1.3.5. Weak solutions to (1) and their local existence. We start this section 
with 


Definition 3. Let /r, u > 0 and 0 < T < oo be a fixed positive parameter. 
The function 

u G L°°([0,r) X LL) 
is a (very weak) solution of (1) if 


f f [dt(p-nA°‘(l)+fj.A^dt(p]u+da;(j) 
Jo Jq 


u 


U‘ 





m 


\t=0 '^0, 


for every test function 4>{x,t) G x fl). 

If a solution u verifies the previous definition for every 0 < T < oo, it is 
called a global solution. 


Lemma 1. Let uq G , s > 2 be the initial data for (1) with //, u > 0 and 
0 < a,/3 < 2. Then there exist T{uo,M) and the unique solution 

ugC{[0,T{uo,M)],H^{Q)) 

to (1). The maximal time of existence T{uq,M) is characterized by 

limsup ||u(t)||^i,oo(Q) = oo. 

t^T{uo,M) 


Proof. Using the same ideas as in [24, Theorem 3.1], we can construct solu¬ 
tions to the regularized problem 




U 


'6,5 


2 ^ + M(1 - 


-vA'^Ue^s-nAl^dtUe^s+edlue^s+ddldtUe^s, 


with initial data rte,5(0) = uq. Standard energy estimates give us uniform 
bounds. Then we can pass to the limits e, J —)• 0. The proof of the continu¬ 
ation criteria can be obtained by energy methods. □ 


2. The entropy inequalities 

In this section we provide the proof of three entropy inequalities that, in 
our opinion, may be of independent interest. 

Proposition 4. Let u be a given function and 0<a<2, 0<e<a/2 be 
two fixed constants. Then 

(22) \\u\\%c/ 2 -e,i <C{a,d,e)\\u\\LiXf[u], 


(23) 


\u\ 


< 


Ha/2 — ^ 


olA 


\u\\L^Zf[u], 


provided that the right hand sides are meaningful. 
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Proof. Let us fix i = 1. First we symmetrize 

^ ^ f f u{x)-u{y) 


x?\u] = 


Cm e / p-v. [ I 

-'T'i JTd\x-y- 27 r 7 |'<+“ 

[ RV. f , 

’ A'i Jjd \x-y- 27 rj\^- 


log{u{x))dydx 


E / A V. 


Ijd \x — y — 27r7|M« 
w(x) - u(y) 


log(«(y))(i2/(ix 


TTd J'j'd \x - y - 27rjf 


> 0 . 


Furthermore, since (a — b) log (|) > 0, every term in the series is positive, 
i.e. for every 7 G Z'^, we have 


[ PY [ - u{y) 

Ijd ■ ' Jjd \x - y - 2Tr'y\<^- 


dydx > 0 . 


In particular 

Let us consider first the case 0 < n G We have 
iL.ii _ f f l«(^) 




Jjd Jjd \x — 2 ^ 

^ f f y Xx) - u{y)\ \x - 
JTdJjdJo !_ |x - |x - 

Isrtlx) + (1 — s)n(y)d/^l , , , 
x -^—^dsdxdy 
|srt(x) + (1 — s)u(y)|^/2j 

^ rO.5 rO.5 

— -'1 -'2 ) 


f f W{x) - ujy)]"^ _1_^ 

hdJjdjQ |x-y|M« \su{x) + (1 - s)u{y)\'^ 

f f y\su{x) + {l-s)u{y)\ 
h = I / / -^- ti-y; - dsdxdy. 


dsdxdy, 


Ijd Jjd Jq 


This latter integral is similar to the Riesz potential. Due to the positivity 
of M, we have 

r f f y suix) + (1 — s)u(y) , , , 

= C{e)\\u\\Li, 


C{e) = 


uTTedy- 


We have 


h = / P.V. 


m(x) - u{y) (u{x) 


'^d^-'-hd |x-yX XX) 


dydx. 
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Consequently, we get 




C 


a,d 


The case 0 < n E was first proved by Bae & Granero-Belinchon [2]. For 
the sake of completeness, we include here a sketch of the proof. Using (24), 
we have 


inu] > 


c. 


a,d 


1 


f f \u{x)-u{y)\^ 
2 2\\u\\l^ Jjd Jjd |x-y|'^+" 


■dxdy, 


so 




< 


— P , 
^a,d 

The proof for the case f = 2 is similar. 


|u||LooXf [u]. 


□ 


3. Proof of Proposition 1: Weak solutions 


We prove the result for U = T, but the same proof can be adapted to deal 
with U = M. We consider the regularized problems 

„,2 


(25) dtUe + dx 


= —+ ed^Ue, X E T, t > 0, 


+ M(1 + e — Ug)^ 
with the regularized initial data 

ng(x,0) = Je* uo(x) + e, 

and e < 1/2. These approximate problems have global classical solution 
due to the Theorem 3.1 in [24]. Consequently, we focus on obtaining the 
appropriate e-uniform bounds. By assumption 0 < u(x,0) < 1. We apply 
the same technique as Cordoba &: Cordoba [11] (see also [7, 13, 8, 2, 1, 16, 
17, 18, 5, 19] for more details and application to other partial differential 
equations), i.e. we track 

\\Ueit)\\L^ = Ue{xt,t) = M,(t) 

and 

minug(x,t) = Ue{x^,t) 

Due to smoothness of Ug, we have that A4g(t), 911g(t) are Lipschitz, and 
consequently almost everywhere differentiable. Hence, using 

■M'eii) = ^lk(i)llL- = = ^(Ug(Xi,t) 

together with the kernel expression for A“, we have the e-uniform bounds 

0 < e < Ug(x, t) < 1.5. 

By space integration of (25) we get 

he{i)\\L^ = lke(0)||Ll < 


Ll +7r. 


We compute 

d ^ r 


/ dtUelog{Ue)dx 

Jt 

I 


UedxU. 


+ M(1 -be — Ug)2 


dx — uX“[ug] — eXf [ug]. 
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We have 


, ^ 2 ^ 2ued^u^{l + M)-2M{l + e)d^u^ 

d^\og[u, + M{1 + e - Ue) ) = - 2 I , - 12 -’ 

^ + m( 1 + e — lie)"' 


and, as a consequence 
UcSxUf 


I 


li^ + M(1 + e — lie)^ 


dx = 


M(l + e) 


dxUf 


M + 1 jf u^ + M{1 + e — lie)^ 
M(1 + e) f dxUe 1 

A ui I 


dx 


M + 1 


-M r 

^ ^ Jt 


dx. 


+ M 


l + e — Ue 


M + 

-\/M 
M + 1 

0 . 


+ M 


l + e-Me 


1 + e —Me 




rdx 




arctan 


fVMfi 


+ e - lie 


V V 


Thus, we conclude 

J'i[iie(t)] + / X"[iie(s)]ds + e / Xi[iie(s)]ds = J’i[ne(0)]. 

JO JO 

As 

-Tl[iXe(0)] < C, 

we get a e-uniform estimate. Now we apply (23) from Proposition 4 and 
Ue{x, t) < 1.5 to get 

/ ||w(s)||^c/2ds < C". 

JO 

Consequently, we have e-uniform bounds 

Me(i) G T°°(0, oo; L°°) n L^{0, oo; 

and the first entropy inequality (4). For the second entropy inequality (5), 
we compute 

d 


dt 


= / dtulog [u^ + M{1 — u)'^) dx 

Jt 

u^dtu -|- uM(l — u){—dtu) 


I 

Jt 


+ 2 


+ M(1 - u)2 


dx 


= h+h 

We are going to handle both integrals separately. We have 


h = 


2udxU + 2M(1 — u){—dxu) 


dx 


f V? 

A yj + M(1 - uY u2 + - m)2 

— V f A^ulog [u^ + M{1 — uY) dx 

Jt 

2 — 12 *^^ “ ^ M(1 - u)^) dx 

+ M{1 — uY Jf ' 


- I 

= —u f A'^u log (u^ + M{1 — uY) dx, 
JT 


dx 
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where in the second term we have used 

2udxU + 2M(1 — u){—dxu) 


{v? + M{1 - uYY 
The second integral reads 


= -dx 


1 


h = - 


2v{M + 1) / 




uu 


v? + M{l-uY 


dx — M 


u2 + M(l-u)V ' 

2udtu 


I 


V? + M(1 - u)2 


dx. 


Now we compute 

A 

dt 
so 

2udtu 


log + M(1 - uf ) dg = 2 /■ 

^ j-f u"' + M(1 — 


:dx = 


9 ,,,, ^ — I log iu^ + Mil — uS^) dx 

■u2 + M(l-u)2 l + M[dtyT ^ 

(1 


I 

L arctan ^ dx 

Collecting all these computations we get 

J^3[u{t)] + U f Xg [M(s)]ds + -^ 

Jo J- 


+ M 


—2\/M / arctan ( \/M f 
JT 

M 


log (u(t)^ + M(1 — u(t))^) dx 

/I 




— 1 dx 


= -^3 [■Wo] + 


1 + M 


Ljt 


\uit) 

log (ug + M(1 — uo)^) dx 


-2Vm [ 

Jt 


arctan — 1) ) dx 


4. Proof of Proposition 2: Decay estimates 

Let us prove first the periodic case. The norm is preserved. Conse¬ 
quently, the mean propagates. We again apply the technique of tracking 
A4(t) and Recall 

\\u{t)\\L^ = u{xt,t) = 

The smoothness needed to proceed with JlA{ty is, for this proposition, an 
assumption. 

Since x,y E T, |x — < (27r)^+“. Hence, using (11) and (13), we get 

, 2r(l + a)cos((l-a)7r/2) , , , , ,, 

A“n(xi) > ^- >_±M^u{xt,t) - (uo)). 


( 271)1 


Consequently, 


d 


M'{t) = — {u{xt,t) - (uo)) 

= -A".R) 

2r(l + a) cos((l — Q;)7r/2) 


< — 


(271) 


l+a 


iu{xt,t) - (uo)). 
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Integrating this ODI, we have 

2r(l + Q:) cos((l —a)7r/2) , 

II«WIIl°°(T) < (uo) + (||wo||l°° - {uo))e . 

Let us turn our attention to the flat at inflnity case. Again the norm 
propagates. We take a positive number r > 0 (that will be specified below) 
and define 

Ui = {r] £ B{0, r) s.t. u{xt) - u(xt - rj) > u{xt, t)/2}, 
and U 2 = B{0, r) —U\. We have 

ui^t) 


Li= \u(xt - y)\dy > / u{xt - y)dy 
JM. J IA 2 

so (for u(xt) > 0), 

2|ko||Li 


> 


■\U2l 


(26) 


u{xt) 

A°‘u{xt) = c(a)P.V. J 

> c(a)P.V.^ 
u{xt) 


< -m. 


u{xt) - u{xt - y) 

\y\l+a 

u{xt) - u(xt - y) 


Ui 


> c{a) 


2r- 


We choose now 


r = 


2||no||Li 
u(xt) 


\y\l+a 

2||uo||li 

u(xt) 


dy 

dy 


thus, recalling (13), we have that 

A°‘u{xt) > 


Cq:,i u(xt)^+" 


2'+“ ll^o|l2i 

for both u{xt) > 0 and u{xt) = 0. With the same argument as in the 
periodic case, we get 

du ||u(t)||itf 

thus, using explicit value of Ca,i we arrive at 


|^^(^)||L°° < 


i + a IholPi V 


5. Proof of Theorem 1: Global solutions for max{a,/3} > 1 

Equipped with the Lemma 1 and its proof, we can focus on the appropriate 
energy estimates that ensures global existence (rigorously, we should do this 
on the level of the regularized problem) 

Notice that we also have a global bound 

(27) \\u{t)\\l2 + y\\u{t)f^p/2 + 2u [ ||u(s)||^„/2ds < ||uo||l 2 +/i||uo||^;3/2- 

Jo 














ON THE GENERALIZED BUCKLEY-LEVERETT EQUATION 


15 


We split the proof in three parts; the first one is devoted to the proof of the 
purely parabolic case /r = 0. Then, in step 2 and 3 we consider the cases 
/i>0/3<l<Q; and /U>0a<l</3, respectively. 

Step 1: Case ^ = 0. We perform the estimates for s = 1, the case s > 1 
being analogous. Testing (1) against Art and using the self-adjointness, we 
have that 

1 d „ 


(28) 

with 

(29) Ii = - 


2^^ 11 ^^ 11 ^ 0.5 


= h + h, 

v?{2u — 2M(1 — u)) 
v? + M{l-uY {v?‘ + M{l-uYf 


2u 


druAudx 


(30) 

Due to Holder’s inequality 


h = 1+0 

H^r- 


(31) 


h < ||it| 


2u 


v?{2u — 2M(1 — u)) 


u‘^ + M{l-uY {v? + M{l-uYf 

Using interpolation 


jjl < c 


i-i 

a 

^0.5 


H- 


1 + 0 ) 


and (8) we have 


Ii < \\u 


|2-- I 

a 

1^0.5 I 


1“ i+oC(rto,M) < \\uYho.5c{uo,M,Y + u||tt|p. i+„ . 


We conclude 

(32) ||n(t)||H0.5 < 

Testing (1) against —d'^u and integrating, we have 


(33) 
with 

(34) h = - 

(35) 


Id 2 


2u 


— h + h, 


m2 + M(1-'u)2 {u‘^ + M{l-uf) 


uH2u-2M(l-u)) \ 
(n2 + M(l-n)2)2 y 


h = -u||u||^i+o. 


We have the bounds 


(36) 


/ 2m 

Y{2u-2M{1-u))\ 

\m2 + M{1 — m)2 

(m2 + M(1 -m)2)2 ) 

/ 2m 

m2(2m — 2M(1 — u))\ 

V'a2 + M(l-M)2 

(m2 + M(1-m)2)2 J 


< c{uo,M)\\u\\l 2 , 


L2 


< c{uo,M)\\u\\hi, 


m 


so, using interpolation, we have 

2u {2u — 2M {1 — u)) 

v? + M{l-uY (m2 + M (1 - rt)2)2 


< c{uq,M)\\u\\^^-c,/2. 


_H'l-o/2 
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We use the duality pairing “/2 _^o/2 i together with Moser’s inequality 
and embeddings to obtain 

^ — ZjIVI — Ubj j \ ^ 

r,U 


h < 

< C 


2u 


+ c 


v? + M{l-uY 

(m2 + M(1 - uff ) 

f 2u 

m^(2m — 2M(1 — m)) \ 

\u^ + M{l-uY 

(m2 + M(1 -m)2)2 ) 

f 2u 

m2(2m-2M(1-m))' 

\v? + M{l-uY 

(m2 + M(1 -m)2)2 , 


\\O^U\\fjo./2-l 

Hl-a/2 


Hl-a/2 


115x111^°° ||'u||j:/1+Q!/2 




Notice that the above estimate can be also derived using the Kato-Ponce 
inequality. As a consequence, 


^ Jo 

Step 2: Case /i > 0, /3 > 1 (and 0 < a < 2). We consider the case 
s = the other cases (s > being analogous. Using (27) together 
with /3 > 1, we have 


1111 (^) 111 .°° < C\\u{t)\\HP/2 < C{l3,n,M,uo). 

Now we test (1) against Am. Due to (29)-(31) with (8), we obtain 

Integrating, we obtain 

h{t)\\% 0.5 +^J-\\uit)f. 1+/3 + [ 2zi||M(s)f y^ds < C\\uof 

Testing against —d‘^u we can perform energy estimates as in Step 1. We 
obtain 

We use the interpolation inequality 

\\f\\L^<c\\f\\f,/2\\Wj:^. 


so 


C\\dxU{t)\\l3 < C'||da;M||^4||dxll||L2 

Now we can use Gronwall’s inequality to obtain 






H 


1+1 


|M(s)f.,^.ds < C(/3,/r,M,Mo)e^(^’^’^’““)*. 

H 
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Step 3: Case ^>0, 0</3<l and u > 0, 1 < a < 2. As before, we test 
(1) against Art. Using (8), we obtain 

Now we use the interpolation 


|n(t)||^i <C||uo||^V||u(t)|| : 

H 2 


2-/3 
I l+a —/5 


For q: > 1 we have that 


hence 


2-/3 
1 + a - /3 


< 1, 


llllolP 1+/9 
H 2 


II^^WII^o.5 + 1+/3 + 2u / ||M(s)f 

Jo 

Now, we test against —d^u. We can conclude as in Step 1. We obtain 


W{t)\\%i + n\\u{t)f^^_^p +u \\u{s)f^,+^ds 




6. Proof of Theorem 2: Global solution for maxla,/3} = 1 
Step 1: Case u > 0, a = 1, // = 0. 

We do the case s = 1, the other cases being analogous. Testing (1) 
against Au and using the self-adjointness, we have equations (28),(29), (30) 
and (31). Notice that, under the hypothesis 

\\uo\\l^ < 7- 

Since ||tt(t)||Loo < ||rto||L°°) we have that 

u^(2u — 2M(1 — u)) 


2u 


u'^ + M{l-uf (ri2+ M( 1 -n)2)2 


^ 27(M + 1) ^ 27^( 7 + M)(M + 1)^ 


Using that 

we conclude 
(37) 


M 


27(M + 1) 27^(7 + M){M + if ^ ^ 


M2 


M 


M2 


|n(t)|| 


HO.5 


+ 6 


f lk( 

Jo 


s)||^ids < 


Ho.s, 


for a small enough 0 < <5. Notice that this 6 only depends on M, uq and u. 

Next, testing (1) against —d'^u and integrating by parts, we get (33), (34) 
and (35). If we integrate by parts in (34), we get 


(38) h < c{uo, M)\\u\\Hi\\dxu\\l4 < c(uo,M,u)||'u||^i -||u|||.i .5 -h c. 

The first inequality above uses also (36) and the second the interpolation 

il< < cii/r„Lii/ii»i 
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Using (38) in (33), we obtain 

^lk( 0 ll^i +z^lk(i)ll^i .5 < c(uo,M,zy)||n||^i +c, 
and, due to Gronwall’s inequality together wit (37), we obtain 

\W\\h^+^ f \\u{s)\\jji.sds < 

Jo 

This ends the proof of case (ii) of our thesis. 

Step 2: Case > 0, a = 1,/r > 0,/3 < 1. In this case we can not use 
the pointwise methods, so we cannot get immediately ||rt(t)|| 1^00 < Ikoll 
Estimate (27) implies a global bound in but this bound is too weak 

to give us a pointwise estimate for u. However, as 

1111011^°° < C'slkoll i+£ < 7) 

H 2 

testing ( 1 ) against Au and using the definition of 7 and 7 *, we have, as in 
step 1 , 

(ll«(«)ll*o,, + ) + «ll«(«)ll*. < 0. 

As a consequence, we obtain the global bound 

\Ht)\\%o.s + f^\\u{t)f i+p + 6 f ||n(s)||^ids < ||iio|||o.6+/^lkof. 1+/3- 
H^r- Jq h^t- 

Now we test against —d'^u and we conclude as in Step 1. Case (iii) is proved. 

Step 3: Case z/ > 0, /i > 0, /3 = 1. In this case, (27) implies a global 
bound in . Then, testing ( 1 ) against Au and using ( 8 ), we have 

\jt ^ C{M)\\u{t)f^,. 

As a consequence, we can apply Gronwall’s inequality to get a global bound 

Mt)\\%o.s + + u j \\u{t)f ^ < C(Ai,M)e‘^(*^)*||uo||^i. 

Jo H^r- 

Now we test against —d'^u and we conclude as in Step 1. Case (i) is proved. 


7. Proof of Theorem 3: Global solution if 0 < a < 1 and /i = 0 


We consider the case s = 2, the other cases being similar. Let us write xt 
for the point where dxU reaches its maximum, i.e. 

dxu{xt,t) = m.gLK{dxu{x,t)) = Mi{t) > 0. 

X 

With a similar argument as in the proof of Proposition 2 (see also [11]), we 
have 

u‘^{2u-2M{l-u)) 


jL 

dt 


Mi{t) = -dx 


2u 


ii2 + M(1-ii) 2 (n2 + M(l-w)2)2 

Due to the kernel expression (11) and (13), we have 
-i>A°'dxu{xt) < -vCa,idxu{xt)- 


drU 


—vA^dxU 


x=xt 


x=xt 
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Due to the smallness choice (10) and ||ix(t)||ioo < ||no||Loo, we have 

2m m^(2m — 2M(1 — m)) ' 


A = 




m2+M(1-m)2 (m2 + M(1-m)2)2 


< S(7*) 

< uCq !• 


Consequently, 






2m 


M' 


v?{2u — 2M(1 — m)) 
+ M(1-m)2 (m2+ M(1-m)2)^ 


— uC, 


0,1 


OtU < 0. 


Let us write for the point where dxU reaches its minimum, i.e. 
dxu{x^,t) = xamdxu{x,t) = 

X 

As before, due to the kernel expression (11) and (13), we have 
— l'A'^dxU{Xf) > —vCa,idxu{xi). 

Consequently, with the same argument, we have (for negative dxu{x^,t)) 


= - 8 , 


2m 


m2(2m — 2M(1 — m)) 


dn 


u 


> -dxu{x-t) 

> 0 . 


m2+M(1-m)2 (m2 + M(1 - m)2)^ 

2m m^(2m — 2M(1 — m)) 


dx 


m2 + M(1-m)2 (m2 + M(1 - m)2)^ 


— vK°'dxu{xi) 
+ uCa 1 


Hence 


(39) ||M(t)|| 

We test equation (1) against 5 ^m and integrate by parts. We have 


(40) + 

with 


(41) ^5 =( 


2m 


\v? + M{1 - uf 


(42) le =-u||m||^2+:|. 


m^(2m — 2M(1 — m))\ 
(m2 + M(1 — m)2)^ / 



dxUdx 


Using (39), we have 

h < ciuo,M)\\uf^2, 

thus, applying Gronwall’s inequality, we have 





















20 


J. BURCZAK, R. GRANERO-BELINCHON, AND G. K. LULI 


8. Proof of Proposition 3: Finite time singularities 
First, we study the case u = 0. Let us take uq such that 


uo > 0, uo(0) = 0, 


and 

(43) J„ = r < oo. 

J— I l®l 

We argue by contradiction: assume that we have u{t) a global solution 
corresponding to uq. Recalling the expression a{x) given in (21) we define 
the characteristic curve y{t), solution to 

(44) y'{t) = a{u{y{t),t)), y{0) = 0. 
and 


v{x,t) = u{x + y{t),t). 


Notice that, due to (1), 


Thus 


A 

dt 


u{y{t),t) = 0. 


u{y{t),t) = uo(y(o)) = 0 = u(o,t). 


Now we have 


dtv{x) = dtu{x + y{t),t) + dxv{x)a{v{0)) 
= dxv{x){a{v{0)) - a{v{x))) 

(45) = -d^ ifivix))) ., 

with / given by (20). For a fixed 0 < <5 < 1, we define 

(j){x) = |x|"‘^l[_i,o] 

and 


(46) 


J{t) = / {v{x,t) — v{0,t))4>{x)dx. 

Jm. 


Notice that if J{t) blows up, due to the inequality 

J{t) = f < ||u(t)||c« = \\uit)\\cs, 

J-I \x\ 


the solution forms a singularity. 
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0=0.25,v=M=0.5 



Figure 1. Evolution in the case a = 0.25, u = M = 0.5, = 0. 


Testing equation (45) against 4>{x), we have 

— / v{x)4>{x)dx = — / dx{f{v{x)))4>{x)dx 


dt 


1-1 


'-1 


f{v{x))dx(pix)dx + 


'-1 


= 6 




1 


'-1 

rO 


X 


1+5 


dx + f{v{-l)) 


125 


= <5 / f{v{x))^{xY 


'-1 


11+5 


dx + /(n(-l)) 


> 


> 


1 + M 
5 / 


{v{x)(j){x)Y dx 


1 + M VV-i 


v{x)4>{x)dx 


where we have used the positivity of the solution and Jensen’s inequality. 
We obtain the ODI 


d_ 

dt 


J{t)> 


1 + M 


Jit? 


and the blow up of J{t) in finite time T* = T*{5,uq,M). 

We have proved the case u = 0, but the proof of the case 0 < u and a = 0 
is analogous and can be easily adapted from here. 


9. Numerical simulations 

In this section we present our numerical simulations suggesting a finite 
time blow up in the case u > 0, 0 < a < 1. To approximate the solution, we 
discretize using the Fast Fourier Transform with N = 2^^ spatial nodes. The 
main advantage of this numerical scheme is that the differential operators are 
multipliers on the Fourier side. Once the spatial part has been discretized, 
we use a Runge-Kutta scheme to advance in the time variable. 
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a=1,v=M=0.5 



Figure 2. Evolution in the case a = 1, u = M = 0.5, // = 0. 


In our simulations, we consider the initial data 
(47) uo(x) = l-e-^" (^1-^y 


and values M = v = 0.5. and /r = 0. Then, we approximate the solution for 

(1) for different values of the parameter 0 < a < 1. In particular, we study 
four cases: 


(1) a = 0.25, 

(2) a = 0.5, 

(3) a = 0.75, 

(4) a = 1. 


Qualitatively, the evolution in the cases a = 0.25, a = 0.5 and a = 0.75 
looks alike. The numerics suggests that a blow up of the derivative appears 
in finite time in the cases a = 0.25, a = 0.5 and a = 0.75 (see Figures 1 and 
3): 

limsup \\dxu{t)\\L^ = oo. 


However, in the case a = 1, the solution seems to exists globally (see Figure 

2 ). 

In Figure 3, we plot the evolution of \\dxu{t)\\L°° ■ This figure shows that 
in the critical case a = 1, the derivative may grow for short time, even if it 
remains globally bounded for large times. 

Next, we add the term ^A^dtu. We consider the same initial data (47) 
and values M = u = ^ = 0.5. Then, we approximate the solution for (1) 
for different values of the parameters 0 < a, /3 < 1. In particular, we study 
three cases: 


(1) a = [3 = 0.25, 

(2) a = 13 = 0.5, 

(3) a = 13 = 0.75. 

Interestingly, we observe (see Figure 4) that even for small values of a and 
/3, in the case with /r > 0, there is not evidence of hnite time singularities. 
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Figure 3. Evolution of \\dxu{t)\\L°°, v = M = 0.5, = 0. 



Figure 4. Evolution of \\dxu{t)\\L°°, v = M = ^ = 0.5. 
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